Statistics of the critical percolation backbone 
with spatial long-range correlations 
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We study the statistics of the backbone cluster between two sites separated by distance r in 
two-dimensional percolation networks subjected to spatial long-range correlations. We find that 



the distribution of backbone mass follows the scaling ansatz, P(Mb) ~ M } 



-(a+l) 



f(M B /M ), where 



f(x) = (a+r/x v ) exp(— x v ) is a cutoff function, and Mo and rj are cutoff parameters. Our results from 
extensive computational simulations indicate that this scaling form is applicable to both correlated 
and uncorrelated cases. We show that the exponent a can be directly related to the fractal dimension 
of the backbone ds, and should therefore depend on the imposed degree of long-range correlations. 
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I. INTRODUCTION 



Percolation is a useful model to study a variety of sys- 
tems in many fields of science displaying structural disor- 
der and statistical self-similarity. In particular, the per- 
colation geometry has been frequently used as a simple 
paradigm to investigate fluid flow in porous media [1-3]. 
In this type of problem, the geometry underlying the sys- 
tem can be very complex and display heterogeneous fea- 
tures over a very wide range of length scales, going from 
centimeters to kilometers [4] . For example, an open ques- 
tion in the modeling process of oil recovery is the effect 
of the connectedness of the porous medium on the dy- 
namical process of fluid displacement. If the pore space 
is so poorly connected as to be considered macroscopi- 
cally heterogeneous, one expects the overall behavior of 
the flowing system to display significant anomalies. In 
this way, it is important to investigate the physics of 
disordered media at a marginal state of connectivity, for 
example, in terms of the geometry of the spanning cluster 
at the percolation threshold [5,6]. 

The most relevant subset of the percolation cluster for 
transport is the conducting backbone. It can be defined 
as the cluster that carries the current when a voltage dif- 
ference is applied between two sites. Thus the backbone 
structure alone determines the conductivity of the whole 
percolation network. Recently, Barthclemy et al. [7] have 
shown that the average mass of the backbone (Mb) con- 
necting two sites in a two-dimensional system of size L 
obeys the scaling form (Mb) ~ L dB G(r / L), where the 
function G(r/L) can be approximated by a power-law 
for any value of r/L. Their results from numerical sim- 
ulations with uncorrelated structures indicate that, for 
the case where r ~ L, the distribution of backbone mass 
is peaked around L dB . When r -C L, the distribution 
follows a power-law behavior. 

The fact that the geometry of real rocks is generally 
random does not necessarily imply that their disordered 
morphology is spatially uncorrelated. In other words, 



the probability for a small region in the system to have 
a given porosity (or permeability) may not be indepen- 
dent of the probability associated with other locations. 
This is certainly true for some types of rocks and geolog- 
ical fields where geometrical and/or transport properties 
can be adequately characterized in terms of their spatial 
correlations [3]. It is under this framework that the cor- 
related percolation model [8-11] represents a more real- 
istic description for the pore space in terms of structure 
and transport. In a previous work [12], the hydrody- 
namic dispersion behavior of percolation porous media 
with spatial correlations has been investigated. More re- 
cently, we studied the displacement dynamics between 
two fluids flowing through correlated percolation clusters 
[13]. It was found that the presence of correlations can 
substantially modify the scaling behavior of the distribu- 
tions of relevant transport properties and, therefore, their 
universality class. In the present paper we extend the 
previous work of Barthelemy et al. [7] to investigate the 
effect of spatial long-range correlations on the statistics 
of the backbone mass connecting two "wells" in a two- 
dimensional percolation geometry. In Sec. II, we present 
the mathematical model to simulate long-range spatial 
correlations. The results and analysis of the numerical 
simulations are shown in Sec. Ill, while the conclusions 
are presented in Sec. IV. 



II. MODEL 



Our model for the porous medium is based on a two- 
dimensional site percolation cluster of size L at criticality 
[5] where correlations among the elementary units of the 
lattice are systematically introduced [10,11]. For a given 
realization of the correlated network, we extract the per- 
colation backbone between two sites A and B separated 
by an Euclidian distance r <C L that belong to the infi- 
nite cluster and are sufficiently far from the edges of the 
system to prevent boundary effects. The correlations are 
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induced by means of the Fourier filtering method, where 
a set of random variables tt(r) is introduced following a 
power-law correlation function of the form 



(u(r)u(r + R)) oc R~ 



[0 < 7 < 2] . (1) 



Here 7 = 2 is the uncorrelated case and 7 « corre- 
sponds to the maximum correlation. The correlated vari- 
ables zt(r) are used to define the occupancy C( r ) 01 the 
sites 



C(r) = efe-«(r)] , 



(2) 



where is the Heavyside function and the parameter (f> 
is chosen to produce a network at the percolation thresh- 
old. Due to computational limitations, we restricted our 
simulations to two values of 7 = 2 and 0.5, correspond- 
ing to uncorrelated and correlated percolation structures, 
respectively. For each value of 7, we performed simula- 
tions for 100 000 network realizations of size LxL, where 
L = 1024, and different values of the "well" distance r 
to compute the distribution P{Mb) and the cumulative 
distribution F{Mb) of backbone mass 



/■OO 

F{M B ) = 

J M B 



P(M)dM . 



(3) 



III. RESULTS 

In Fig. 1 we show the log-log plot of typical cumulative 
distributions of backbone mass F(Mb) for uncorrelated 
as well as correlated morphologies. It is clear from this 
figure that F{Mb) displays power-law behavior for in- 
termediate mass values in both cases. In addition, the 
scaling region is followed by a sudden cutoff that decays 
faster than exponential. A similar behavior has been ob- 
served experimentally [14] and through numerical simu- 
lations [15-18] for the phenomenon of impact fragmen- 
tation. Based on these features, we argue that F(Mb) 
should obey the following scaling ansatz: 



F(Mb) ~ Mg a exp 



Mo) . 



(4) 



where a is a scaling exponent and M and 77 are cutoff 
parameters. For comparison with the approach presented 
in Ref. [7] , here we determine these parameters from the 
simulation results in terms of their corresponding dis- 
tributions of backbone mass P{M B )- From Eq. (4), it 
follows that 



P(Mb) ~ MB- {a+1) f{M B /M Q ) , 
where the function /(Mb/Mo) has the form 

'^VlexpL^V 



f {M B /Mo) = 



oe + rj 



M J 



\Mo 



(5) 



• (6) 



Figure 2 shows the distributions P(Mb) generated for the 
uncorrelated case (7 = 2) and different values of r. They 
all display a lower cutoff of order r (the smallest back- 
bone possible is a straight line connecting points A and 
B) and an upper cutoff of order L Db , where a "bump" 
can also be observed [7]. For comparison, each of these 
distributions has been rescaled by its corresponding value 
at the position of this bump. The smaller the distance 
r is between the wells, larger is the range over which 
the scaling term of Eq. (5) holds, P(M B ) ~ M B ~ (a+1) ■ 
The solid line in Fig. 2 corresponds to the best nonlin- 
ear fit we found for the r = 8 data in both the scaling 
and cutoff zones with a = 0.255, M = 3.32 x 10 4 and 
r\ = 1.5. Because the sampling of the backbone mass 
should be equivalent to the sampling of blobs in the per- 
colation cluster [19] , it is possible to draw a direct relation 
between the scaling exponent a and the fractal dimen- 
sion of the backbone [20] . Accordingly, the exponent r 
governing the blob size distribution can be calculated as 
t = d/e?B + 1. Since the exponent tb = a + 1 governs the 
statistics of the entire backbone, we obtain by integra- 
tion that a = d/ds — 1, which gives a fractal dimension 
of cIb ~ 1.6. This is in good agreement with the current 
numerical estimate of d B = 1.6432 ± 0.0008 [21]. 

We now turn to the case with spatial long-range corre- 
lations. As shown in Fig. 3, Eq. (5) (solid line) also fits 
very well the simulation data for P(Mb) with 7 = 0.5 
and r = 16, both in the scaling region and in the cutoff 
zone. The parameter set obtained with a nonlinear es- 
timation algorithm includes a = 0.075, M — 5.4 x 10 5 , 
and 7] = 2. Compared to the uncorrelated geometry, the 
results shown in Fig. 3 for large masses clearly indicate 
the presence of a narrow plateau followed by a much more 
pronounced bump and a sharper cutoff profile. These 
features are consistent with the exponent a ss and also 
reflect the fact that 77 is significantly larger for the corre- 
lated case. The resulting fractal dimension of the back- 
bone dB ~ 1-86 is in good agreement with previous es- 
timates for correlated structures generated with 7 = 0.5 
[10,13]. 

For completeness, we show that the scaling ansatz used 
to fit the simulation results for P(Mb) is consistent with 
the observed behavior for F(Mb) over the whole range 
of relevant backbone masses. Precisely, the solid lines in 
Fig. 1 correspond to Eq. (3) with the same set of param- 
eters used to fit the simulation data for 7 = 0.5 and 2 
in terms of Eq. (5). Although close to zero, the value of 
a = 0.075 for the correlated geometry is sufficiently large 
to characterize the power-law signature of the cumulative 
mass distribution F(Mb)- The differences in the expo- 
nent a obtained for correlated and uncorrelated cases can 
be explained in terms of the morphology of the conduct- 
ing backbone. As 7 decreases, the backbone becomes 
gradually more compact [10]. This distinctive aspect of 
the correlated geometry tends to increase the value of ds 
and therefore reduce the value of a as the strength of the 
long-range correlations increases (i.e., 7 decreases). As 



2 



in Rcf. [7], here we also investigate the scaling behavior 
of the average backbone mass (Mb), but for a correlated 
geometry. In Fig. 4 we show the log-log plot of (Mb) 
against r for 7 = 0.5 and three different values of L. 
The approximation proposed in Rcf. [7] for uncorrelated 
networks when r«i reveals that (Mb) should scale as 

(M) ~ L dB -^r^ , (7) 

where the exponent ip is the codimension of the fractal 
backbone, i.e. ip = d — d B - Using Eq. (7) and d B = 1-85, 
the inset of Figs. 4 shows the data collapse obtained 
by rescaling Mb and r to the corresponding values of 
L ds and L, respectively. From the least-squares fit to 
the data in the scaling region, we obtain the exponent 
ip = 0.17 ± 0.03. This result is in good agreement with 
the estimated value 0.15 for the fractal codimension. 

Finally, it is also interesting to investigate the case 
where r « L. As for uncorrelated clusters [7], we expect 
the distribution P(M B ) for correlated geometries to obey 
the simple scaling form 

1 / Mb \ 
P(M B ) ~ j^g [j^j , (8) 

where g is a scaling function. In Fig. 5 we show re- 
sults from simulations with correlated networks gener- 
ated with 7 — 0.5 for r « L and three different values 
of L. It can be seen that the probability distribution 
is peaked around an average value (Mb) of the order of 
L dB . As depicted, Eq. (8) can indeed be used to collapse 
the data. The value adopted for the fractal dimension of 
the backbone, ds = 1-86, is compatible with the corre- 
sponding degree of correlation imposed. 



IV. SUMMARY 

We have studied the scaling characteristics of the 
backbone mass distribution between two sites in two- 
dimensional percolation lattices subjected to long-range 
correlations. A scaling ansatz that is capable of describ- 
ing the power-law region as well as the complex details 
of the cutoff profile is proposed and it is shown to be ap- 
plicable for both correlated and uncorrelated structures. 
Based on the results of extensive simulations, we find 
that the presence of long-range correlations can substan- 
tially modify the scaling exponents of these distributions 
and, therefore, their universality class. We explain this 
change of behavior in terms of the morphological dif- 
ferences among uncorrelated and correlated pore spaces 
generated at criticality. Compared to the uncorrelated 
case, the backbone clusters with spatial long-range cor- 
relations have a more compact geometry. The level of 
compactness depends, of course, on the degree of corre- 
lations 7 introduced during the generation process. 



V. ACKNOWLEDGMENTS 

This work has been supported by CNPq, CAPES and 
FUNCAP. 



[1] J. Bear, Dynamics of Fluids in Porous Materials (Else- 
vier, New York, 1972). 

[2] F. A. Dullien, Porous Media - Fluid Transport and Pore 
Structure (Academic, New York, 1979). 

[3] M. Sahimi, Flow and Transport in Porous Media and 
Fractured Rock (VCH, Boston, 1995) and the extensive 
references therein. 

[4] P. R. King, in North Sea Oil and Gas Reservoirs III, 
edited by A. T. Buller et al. (Graham and Trotman, Lon- 
don, 1990). 

[5] D. Stauffer and A. Aharony, Introduction to Percolation 

Theory (Taylor Francis, Philadelphia, 1994). 
[6] Fractals and Disordered Systems 2nd ed., edited by A. 

Bunde and S. Havlin (Springer- Verlag, New York, 1996). 
[7] M. Barthelemy, S. V. Buldyrev, S. Havlin, and H. E. 

Stanley, Phys. Rev. E. bf 60, R1123 (1999) 
[8] S. Havlin, R. B. Selinger, M. Schwartz, H. E. Stanley, 

and A. Bunde, Phys. Rev. Lett. 61, 1438 (1988). 
[9] S. Havlin, M. Schwartz, R. B. Selinger, A. Bunde, and 
H. E. Stanley, Phys. Rev. A 40, 1717 (1989). 
[10] S. Prakash, S. Havlin, M. Schwartz, and H. E. Stanley 

Phys. Rev. A. 46, R1724 (1992). 
[11] H. A. Makse, S. Havlin, M. Schwartz, and H. E. Stanley 

Phys. Rev. E. 53, 5445 (1996). 
[12] H. A. Makse, J. S. Andrade, Jr., and H. E. Stanley, Phys. 

Rev. E. 61, 583 (2000). 
[13] A. D. Araiijo, A. A. Moreira, H. A. Makse, H. E. Stanley, 
and J. S. Andrade Jr., Phys. Rev. E 66, 046304 (2002). 
[14] A. Meibom, and I. Balslev, Phys. Rev. Lett. 76, 2492 
(1996). 

[15] H. Inaoka, E. Toyosawa, and H. Takayasu, Phys. Rev. 

Lett. 78, 3455 (1997). 
[16] A. Diehl, H. A. Carmona, L. E. Araripe, J. S. Andrade 

Jr., and G. A. Farias, Phys. Rev. E 62, 4742 (2000). 
[17] J. A. Astrom, R. P. Linna, and J. Timonen, Phys. Rev. 

E 65, 048101 (2002). 
[18] A. Diehl, J. S. Andrade Jr., and G. A. Farias, G. A., 

Phys. Rev. E 65, 048102 (2002). 
[19] H. E. Stanley, J. Phys. A: Math. Gen. 10, L211 (1977). 
[20] H. J. Herrmann, and H. E. Stanley, Phys. Rev. Lett. 53, 

1121 (1984). 
[21] P. Grassberger, Physica A 262, 251 (1999). 



3 
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FIG. 1. Log-log plot of the cumulative distribution of back- 
bone mass for uncorrelated (7 = 2, circles), and correlated 
networks (7 = 0.5, squares). The solid lines correspond to 
the scaling function Eq. (4) with the same values of the pa- 
rameters obtained from the best fit to the data of Eq. (5) . 
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FIG. 3. Log-log plot of the distribution of backbone mass 
of correlated percolation networks for L = 1024 and r = 16 
(circles), 32 (triangles up), 64 (squares), and 128 (trian- 
gles down). The solid line is the best fit of Eq. (5) to the 
data for the scaling region and lower cutoff, with parameters 
a = 0.075, r) = 2.0, and M = 5.4 x 10 5 . The degree of 
correlation is 7 = 0.5. 
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FIG. 2. Log-log plot of the distribution of backbone mass 
of uncorrelated percolation networks for L — 1024 and r — 8 
(circles), 16 (triangles up), 32 (squares), 64 (triangles down), 
and 128 (diamonds). The solid line is the best fit of Eq. (5) 
to the data for the scaling region and lower cutoff, with pa- 
rameters a = 0.255, Tj = 1.5, and M = 3.32 x 10 4 . 




FIG. 4. Log-log plot of the average backbone mass (Mb) 
for correlated structures (7 = 0.5 against the distance r for 
L = (circles) , (triangles) and (squares) . The inset shows the 
best data collapse obtained by rescaling Mb and r to L 1 ' 85 
and L, respectively. The least-squares fit to the data gives 
the scaling exponent tp = 0.17 ± 0.03. 
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FIG. 5. Data collapse of the distribution of backbone mass 
for correlated networks (7 = 0.5) and three values of r ~ L. 
For each value of L, 100 000 realizations have been generated 
to produce the statistics. Equation (8) with ds = 1.86 has 
been used to collapse the results. 
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